f{R) gravity from the renormalisation group 
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We explore the cosmological dynamics of an effective f{R) model constructed from a renormal- 
isation group (RG) improvement of the Einstein-Hilbert action, using the non-perturbative beta 
functions of the exact renormalisation group equation. The resulting f{R) model has some remark- 
able properties. It naturally exhibits an unstable de Sitter era in the ultraviolet (UV), dynamically 
connected to a stable de Sitter era in the IR, via a period of radiation and matter domination, 
thereby describing a non-singular universe. We find that the UV de Sitter point is one of an in- 
finite set, which make the UV RG fixed point inaccessible to classical cosmological evolution. In 
the vicinity of the fixed point, the model behaves as R^ gravity, while it correctly recovers General 
Relativity at solar system scales. In this simplified model, the fluctuations are too large to be the 
observed ones, and more ingredients in the action are needed. 

PACS numbers: 04.50.Kd,98.80.-k 



I. INTRODUCTION 

Einstein gravity, although successful at solar and galac- 
tic scales, is chaUenged by cosmological observations of 
the early and late time Universe. Both early (inflation) 
and late time acceleration of our Universe require either 
the introduction of an extra degree of freedom in the ac- 
tion, like a scalar or tensor field, or a modification of 
gravity itself. In particular, modifications of the gravita- 
tional action include scalar-tensor actions or non-linear 
extensions of Einstein-Hilbert action, probably the most 
famous being the so called Brans-Dicke [1] and f{R) the- 
ories respectively flHZ]- A common characteristic among 
scalar-tensor and modified theories alike is that they lead 
to a modification of Newton's constant Gn, which ac- 
quires a scale dependence, for example in scalar-tensor 
theories through the coupling of gravity with a scalar 
field. 

At the same time, GR with a cosmological constant A 
has been very successful in describing the late time ac- 
celeration of the Universe from a phenomenological point 
of view, but it is unable to account for a primordial in- 
flationary era. Probably, the most challenging problems 
a cosmological constant faces from a theoretical point of 
view is the order of magnitude problem, i.e why it has 
such a tiny value, as well as the coincidence problem, or in 
other words why it is only at recent times that A becomes 
dynamically relevant. In the context of scalar field or 
modified gravity models, the vacuum energy is replaced 
by a dynamically evolving, effective energy-momentum 
tensor, but this only partly solves the problem, as any 
effective energy-momentum tensor has to reproduce the 
tiny value of A today. 
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In the last years, another modification of Einstein grav- 
ity has been suggested in the context of the Asymptotic 
Safety scenario |8j. According to this, Einstein gravity 
is quantisable as a metric theory in a non-perturbative 
way, provided a non-trivial fixed point exists at high ener- 
gies, under the renormalisation group (RG) [SHTB]. The 
existence of an ultraviolet (UV) fixed point ensures a 
well-defined UV behavior in the standard quantisation 
of gravity. Different investigations have shown that a 
non-trivial UV fixed point exists for the Einstein-Hilbert 
and higher truncations |14H18j . providing strong evi- 
dence that quantum gravity is renormalisable in a non- 
perturbative way. 

The starting point for this particular implementation 
of quantum gravity, is the scale dependent effective aver- 
age action r^ [5^1/] , where k is the cut-off energy, or coarse 
graining scale, which defines the scale above which all 
modes are integrated out, i.e. the functional integral is 
performed on modes whose momentum satisfies p^ <^ k^ . 
The effective action satisfies a functional renormalisation 
group equation [Hj (see also [50]), which describes the 
evolution of Tk[g^i,] as a function of the scale k, giving 
rise to an RG flow on the space of different effective ac- 
tions. 

Although any curvature invariant can be in principle 
included in the effective action, the actual study of the 
RG equations requires that a particular truncation is con- 
sidered. In this paper, we will be focusing in the so- 
called Einstein-Hilbert truncation, where the only cou- 
plings are Newton's G = Gk{k) and the cosmological 
constant A = A/j(fc), with k the renormalisation group 
cut-off scale. The running of G and A changes the cosmo- 
logical dynamics resulting from the action, and has also 
been suggested as a possible resolution to the coincidence 
problem [3TH55]. In particular, Refs P51 - I5^ have studied 
the cosmological consequences of a running Newton's G, 
while Refs [33H35] studied the case of a running cosmo- 
logical constant. Comparison with cosmological observa- 



tions, including supernovae data, has been carried out in 

Refs [MlEMaH]. 

In a cosmological context, it is attractive to think of 
the cut-ofF energy scale k as dynamically evolving with 
cosmic time [2T1 [351 EH] • There are different ways to un- 
derstand this connection in an expanding Universe. Since 
in the effective action modes with momentum p^ 3> k^ 
are integrated out, k defines the energy scale of the the- 
ory, i.e. the typical scale at which the couplings in the 
effective action are evaluated. 

The typical energy of particles in an expanding Uni- 
verse with temperature T at a particular time, is directly 
linked with the expansion; the Universe starts off from 
a hot state and cools down as it expands with cosmic 
time t and in particular, for a homogeneous and isotropic 
Universe described by the Friedman-Robertson- Walker 
(FRW) metric, characterised by the scale factor a{t), the 
typical energy of relativistic particles scales as l/a{t). 
One could then think of identifying the cut-ofF scale k as 



k - ksTit) 



aity 



(1) 



where kg is the Boltzmann constant and Eq a constant 
with dimensions of energy. 

Alternatively, one can think that the horizon size of 
the Universe dH ^ l/H{t), with H = a{t)/a{t) the Hub- 
ble parameter, defines the typical scale of correlations 
between different quantum degrees of freedom, and iden- 
tify 



dH{t)^H-\t). 



(2) 



Notice that both identifications (IT]) and (pi) arc monoton- 
ically decreasing functions of cosmic time, in the context 
of a Hot Big Bang scenario. In addition, both of above 
identifications are performed at the level of the equa- 
tions of motion, for which k{t) is constrained through 
the Bianchi identities pTl HOl - HS] , which provide the con- 
dition for all consistent identifications. 

In this paper, we will consider a different identification 
which is performed at the level of the action, motivated 
by an analogous procedure which generates the effective 
potential for a scalar field theory. We will associate k 
with the scalar curvature, i.e 



e 



R, 



(3) 



through which we can view the RG improved Einstein- 
Hilbert action as an effective f{R) model. The idea 
of Newton's G running with curvature has been sug- 
gested previously HH, although not with the correct 
beta-function, and here we include the cosmological con- 
stant with the full non-perturbative beta functions for 
both couplings. The resulting f{R) model does not in- 
clude the renormalisation effects of matter, or any grav- 
itational invariants other than i?, and so it should be 
viewed as a prototype. However, it will turn out that it 
has some remarkable properties, also allowing us to study 
the RG improved action in an elegant way. One feature 



is that the scale identification is performed at the level 
of the action, in a covariant fashion, so there is no need 
to add extra dynamical conditions through the Bianchi 
identities as described above. 

In the following, we will be working in a unit system 
with c = h = 1, unless otherwise stated, as well as use 
G = rnZ,^ , = Sttk^. Unless otherwise stated, mass scales 
will be presented in Planck units. 



II. RG IMPROVED EINSTEIN-HILBERT 
ACTION 

Our starting point is the RG improved effective action 
in the Einstein-Hilbert truncation. 



rfc[5>] 



f R{9) 



2Ai 



Cr, 



M^g) 



(4) 

with k is the renormalisation group cut-off scale, which 
sets the momentum scale above which modes are inte- 
grated out. The effective, "coarse-grained" action func- 
tional, r/i;[g, ■(/;], interpolates between the true effective 
action in the infrared (IR, fc — >■ 0) and the bare action 
defined in the UV at a cut-off scale /cmax- The inter- 
polation of the effective action as a function of scale is 
controlled by the exact renormalisation group equation 
(ERGE) [19]. If fcmax can be taken to infinity the theory 
is renormalisable, signaled by a UV fixed point in the 
couplings of the theory. 

The quantum corrections can be encoded in the evo- 
lution of the coupling constants as a function of energy,^ 
whose beta-functions can be extracted from the ERGE. 
The form of the latter depends on the choice of the cut- 
off function choice and the gauge. We will follow the 
conventions of |3S1 ST], noting that different choices of 
cut-off function and gauge do not change the qualitative 
features of the beta functions. 

In the standard approach, one defines the dimension- 
less Newton's and cosmological constant as 



g{k) = k^G{k)/2ATT, X{k) = A/k^ 



(5) 



and the running of the dimensionless couplings in d = 4 
is described through the set of first order, coupled differ- 
ential equations El] , 

;^ \ «^ \^- 0\ 10 24g(3s+i(l-3A)) 

dtX = l3x{g,X)^ -2X-12g 2g-|(i-2A)2 ' (6) 

dtg = l3g{g,X)^2g+^^^^^, (7) 

where t = Infc, and l3\, Pg the beta functions. In above 
equations the factor of 247r is included to remove phase 
space factors. 



^ For a cautionary note see I45| . 



There are two fixed points of the above RG flows, a free 
or Gaussian one, with (g* , A*)gfp = (0, 0), and an inter- 
acting one which is attractive in the UV {k — > oo), with 
(5*,A*)uv = (0.015625,0.25) (Refs [g]-[T3 and references 
therein). The existence of a UV fixed point points to con- 
sistent quantum behavior of the system at high energies, 
realising Weinberg's Asymptotic Safety scenario |S]. The 
Gaussian fixed point {k — >■ 0) describes a free theory. 

A phenomenologically viable solution (trajectory) of 
the system ([6])- (I?]) on the <? — A plane is one that starts at 
high energies from the UV fixed point and then evolves 
towards smaller values oi g as k is lowered, passes close 
to the GFP, until it turns to the right towards increas- 
ing values of A. A trajectory passing sufficiently close to 
the GFP will subsequently have a long classical regime, 
i.e G ~ Go, A ~ Aq, with "0" here denoting the present 
value. The classical regime covering many orders of mag- 
nitude in scales is required by terrestrial, solar and galac- 
tic tests, as well as consistency with cosmological evolu- 
tion since Big Bang Nucleosynthesis. 

The Einstein-Hilbert truncation has a couple of fea- 
tures which may not be present in all truncations. The 
eigenvalues of the linearised flow in the vicinity of the UV 
flxed point are complex conjugate, causing oscillatory be- 
havi or of the trajectory around it (see Fig. [I]). Also, the 
flow (6l7) has a singularity at A = 1/2, which terminates 
the classical regime. It has been conjectured that this is 
an artifact of the truncation, and that there may actually 
be a non-trivial fixed point in the IR |39j . 



III. SCALE IDENTIFICATION & EFFECTIVE 
f{R) ACTION 

As explained before, the first step in studying the cos- 
mology of an RG improved action is to identify the cut-off 
scale k with as a function of cosmic time k = k(t). 

In this paper, we will work at the level of the action 
and use a particular ansatz that will allow us to view 
the effects of RG running of the couplings as an effective 
f{R) model, by identifying 



pi?. 



(8) 



where R is the Ricci scalar and p is a dimensionless 
constant. Here, with the particular identification ([8]) 
the dimensionless couplings are defined as g{k) = pR x 
G{R), X{R) = A{R)/{pR), and action Q takes the 
following form 



S 



fiR) 



j4 / 

a xy'—g 
a xy'—g 



R^h{R) 
3847r2 

f{R) , 
2k2 



-h5,„(^,.g) (9) 

5^(^,5), (10) 



with h{R) = p{ 



= ^^l-2pAN 



and the extra factor of 247r ap- 



formed in the beta functions (l6| and (It]). We absorb it 
into the factor k^ = 192'k'^. 

The quantum corrections are now expressed in the non- 
linear effective action, which takes the form of an f{R) 
model ( 10 ). This provides us with a different view of the 
RG effects on the Einstein-Hilbert action Q. What is 
more, the particular scale identification preserves general 
covariance of the action. 

We can compare this procedure with the RG- 
improvement of the effective potential in scalar field the- 
ory [35]. There, if one starts with the tree potential 
V — A(/)^/4!, solves the RG equation for the couphng, 
and makes the identification k = acj), one obtains at one 
loop 



V 



^ ^0 ,4 

A\l-b{Xo)\n{a<l>/koy ' 



(11) 



where 6(A) = SA/IGtt^, which recovers and improves on 
the one-loop effective potential calculated by the stan- 
dard graphical methods. The constants are constrained 
by a renormalisation condition such as 






= K, 



(12) 



where Ar is the physical coupling as inferred, say, from 
a scattering experiment with the background field set at 
(j)o. Normally, we can avoid all mention of a by writing 
V = A0V4!, with A = Ar/(1 - b{X,) ln((/)/0o))- However, 
it is still implicit in the relationship between ^q and the 
scale ko- 

The renormalisation conditions for the effective 
Einstein-Hilbert action can be taken as 



/filflo - 



RfB.-f 



SttGo 



2/fl 



= Ao, 



(13) 



Ra 



where Rq is the curvature scalar evaluated today. 

Some remarks regarding the action ( 10 ) are in order. 



Firstly, we can see that on a fixed point, where h(R) is 
constant, the Lagragian is effectively R^, which is renor- 
malisable [50 . Secondly, there is a singularity of the RG 
flow in the Einstein-Hilbert truncation ^ HD] : the beta 
functions diverge ior Ag = {X — ^)^ . We will therefore re- 
strict ourselves to cosmological evolution which does not 
reach the singularity. 

Finally, let us comment on the dimensionless parame- 
ter p, defined through our identification ([s]), relating the 
RG scale k and the cosmological scale R. We will see in 
Eq. (35 1 that it determines the scalaron mass, and so in 



pearing in the first line because of the rescaling of g per- 



principle could also be fixed. However, as we do not know 
the scalaron mass, we will leave p free, and investigate 
what range of values give an acceptable cosmology. As p 
describes to what extent the RG scale k follows the cur- 
vature R, we would hope to find that p ^ 1: it is natural 
to think of the RG scale as the scale of the important 
dynamics, which in the cosmological context is given by 
the curvature. It will in fact turn out p ^ 1 gives a viable 
cosmology. 




FIG. 1: A viable RG trajectory (blue) on the <? — A plane in the Einstein-Hilbert truncation for the choice p = 1 and the 
initial conditions (37 1. It spirals around the UV RG fixed point and evolves towards the IR as curvature R decreases. The 



intersection of the phase curve with the de Sitter line (black) corresponds to a de Sitter potnt in the cosmological evolution, while 
intersection with the dashed (black) one is where for the slow roll parameter ev ~ 1- The regions where mffj > (rUgfj < 0) 
are seperated by the dotted lines, with m^g the Jordan frame mass squared, defined in (16). The dotted curve consists of two 



separate curves (green and red dots) corresponding to the vanishing of the numerator (denominator) of m^g. They join at the 
upper part of the dotted "ellipsis" , where m^ff remains finite and non-zero. Along the lower part of the "ellipsis" (green) m^ff 
vanishes. The dotted curves outside the "ellipsis" (red) correspond to m^g — >■ oo. Notice that beyond the red dot at A* ~ 0.27 
on the A-axis, m^g becomes negative, and therefore de Sitter space unstable too. 



IV. STABILITY AND THE GR LIMIT 



A. Degrees of freedom and stability conditions 



As a first step to understand the resulting effective 



f{R) action from the renormalisation group (10 1, we 



want to study its stability and its approach to the limit- 
ing case of GR. Below we will introduce some basic facts 
about f{R) gravity that will be necessary for the rest of 
the paper [2H7]. 

It is well known that f{R) models, exhibit an extra, 
massive scalar degree of freedom, dubbed "scalaron" . It 
satisfies a Klein-Gordon type equation, which can be 



found by varying action ( 10 ) with respect to the met 



ric and then taking the trace. 



afniR) 



dVeffjR) 

dfR 



^^T, 



(m), 



(14) 



where D is the d'Alembertian associated to the metric 
5q/3: 2^(111) = 9^'^T(^m)fj.i' is the trace of any matter sources 
present, and 



dV,siR) _ 1 
dfB. ~ 3 



[RMR) - 2/(i?)] , 



(15) 



From ( 15 ) we can deduce the scalaron's mass in the frame 



defined by action ( 10 ) 



_d^V,siR) fR~RL 



RR 



'eff 



dfl 



^fRR 



(16) 



Expression \1Q\ also appears as the effective mass in a 
stability analysis around de Sitter spacetime |5J (and ref- 



erences therein).^ Therefore, stability of the scalaron 
propagation (i.e avoidance of tachyonic instability), as 
well as stability of de Sitter spacetime requires that 



Ks > 0- 



While an unstable scalaron just means that long- 
wavelength scalar fluctuations will grow, the graviton ki- 
netic term must certainly have the correct sign, in order 
to avoid ghosts. This means that 



fR>0, 



(17) 



which at small values of the couplings A and g, is ensured 
through condition 



pR 



>0. 



(18) 



In order to make the connection with the RG, we will 
express both stability conditions, for the scalar and for 
the graviton, in terms of the beta functions, using the 
explicit form of action ^ . 

In the RG-improved Einstein-Hilbert action deriva- 
tives of / can be expressed in terms of RG data, as 



_d_ 
dR 



d_ 
dR 



1 
2R 



4 



P> 



d_ 
dX 



For example, for /r we have 

1 



Ir = 2R 



45 



{hPg + 2p'(ix) 



(19) 



(20) 



while the second derivative is 



fRR = 2h 



45 



-6/i + 2h^ + 2p 



2 Pa 
9 



hPg^g ~ 2p'Px,g 



4.9 



9 



(21) 



Plugging above relatio ns i nto expression (16), and us- 



ing the beta lunctions (6l7l, the scalar mass rr 



'off 



re-expressed as to^jj = ra^{R, g, A). The same is in prin- 
ciple true for fa and other quantities of interest as we 
will see later. 

In particular, from the explicit expression of //? in 
terms of the couplings, one can check that the no-ghost 
condition ( |17[ ) is always satisfied in the domain of inter- 
est, < A < 0.5, < g < 0.02 and of course R> 0. (see 
also Fig. [2] for a plot of fn and Fig. [l] for the phase space 
of a viable RG trajectory.) 



B. The f{R) model in the perturbative regime 

Let us now see how GR is recovered in this framework. 
If an RG trajectory is to be viable, it should have a suffi- 
ciently long classical regime, where any quantum correc- 
tions are suppressed enough not to be observed in astro- 
physical or solar system tests, and therefore the coupling 
constants should be effectively constant, and acquire the 
values observed at these scales. Therefore, GR is re- 
covered in the sense that under the RG flow Newton's 
constant acquires its classical value, G ~ Go = l/"Zp for 
a sufBcient "RG time", large enough to cover the range 
of classical scales (earth, solar and galactic). In the clas- 
sical regime, A has to have a negligible variation too. It 
has been shown that both requirements are achieved if 
the viable RG trajectory passes sufficiently close to the 
GFP at (.g, A) = (0,0) [37]. It is after the close passage 
to the GFP when the classical regime starts, and it turns 
out that the closer the trajectory passes to it, the longer 
it lasts in RG time, and the greater range of scales the 
classical regime covers. 

We will therefore need to first linearise the system of 
beta functions ^-^ around the GFP [27]. To make the 
analysis more clear, it would be better to first proceed 
with the linearisation of the equations without assuming 
any identification for k, i.e keeping k as the independent 
variable in ^-^- We get. 



dtX = -2A + 2ag, 
dt9 = 2ff, 



(22) 
(23) 



with renormalisation group time t = ln(fc/fco), and fco 
a reference scale. The parameter a cut-off function de- 
pendent, but is always positive and of order 1. For the 
optimised cut-off, used to derive the beta functions ^- 
0, a = 6. 



The solution of the linearised system reads 



9 



Krr-i 



1 



A = -aci k 
2 



El ^ Ixt ( — 






(24) 



(25) 



with kT the value of the cut-off scale around the turning 
point in the vicinity of the GFP, and qt = 9{kT), Xt = 
gikx)- Notice also that is. 



Xr/gi 



a, 



(26) 



which implies that At ~ gr, since a ^ 0(1). From 
the above linearised relations we get for the dimensionful 
couplings, 



SttG gT 

— ^^ = ci = ^ = const. 



A — -ac-i k 
2 



C2 = - A 



Tp-+ 2>^TkT, 



(27) 
(28) 



following the notation of Ref. [27]. Equation (27) tells 



us that in this regime, Newton's G becomes a constant, 
and we identify ci = SttGq/k?. For scales k <^ kx, A is 
also effectively constant, and we may identify C2 = Aq. 
Hence 



grXr SttGqAo 



and 



9T 



WttGqAo 



f k^Aq 



From the observed values of Aq and Go, we have 
grj,^ Xj,^ 10-60, kj. - lO^^Vp. 



(29) 



(30) 



(31) 



Let us now turn to the solution of the system under 
the identification k^ = pR. The linearised equations 
are not enough as higher-order terms contribute already 
at 0{R^). An efficient way to include the higher-order 
terms is to substitute into the Talyor expansion around 
the renormalisation point Rq ( 13 ) 



1 



f{R) = /(i?o) + fnln^ (R - Ro) + ^ JrrIr, [R ~ Ro) 



We find that for the optimised cut-off where a = 6, 



(32) 



fiR) 



Go 



[R - 2Ao) + 6(2 - p)piR - Ro) 



(33) 



From the small coupling expansion of ( 16 ) the scalaron 



mass squared in the classical regime is given by 



"^eff,0 



Jr — RJrr. 



3/ 



RR 



Rf) 



Rq 



36(2 -p) g' 



(34) 



(see also Eq. 51 1 and we see that it is positive provided 
< p < 2. Using the renormalisation condition (13) we 
find 



1 



'cff.O 



36(2 - p) SttGo ' 



(35) 



and observe that the scalaron mass is safely at the Planck 
scale, so large deviations from GR at laboratory, solar 
and astrophysical scales are avoided. 

To get an idea of the realistic values of the couplings 
in the classical regime we can evaluate them at solar and 



galactic scales, taking k '-^ R. With R, 



-1/2 



lAU 3 and 



R] 



-1/2 



lO^^m we find that 



5sol 



Rsn] X G. 



sol 



10" 



ffgal 

:i 10-™ 



-Rgal X Ggal ~ 10 , 

(36) 
assuming that Ggai — Gsoi — 10 ^" m^. We see that the 
classical value of the dimensionless coupling g acquires 
a tiny value. For A we cannot follow the same analysis, 
since A has been only measured at cosmological scales, 
k ^ Hq, with Hq the Hubble parameter today. However, 
the product gX ^ GqAq, so Agoi ^ 1. Therefore, the 
values of both g and A on solar and galactic scales lie 
extremely close to the GFP. However, it is intriguing to 
note that by this reasoning, A evaluated at the Hubble 
scale is of order 1, where non-perturbative effects in the 
beta functions are important |27j . 

The form of a phenomenologically viable RG evolution 
on the g — X plane is given in Fig.[l] for the choice of p = 1 
and ^ 



R, 



!n.i„ = 8 X 10-^ i?T = 5 X 10-3, i?,„ax = 50, 
A(i?T) - 10-^ giRx) = 10-^ (37) 



where Rt is the curvature at the turning point close to 
the GFP, in Planck units. The above initial conditions 
are not realistic, but they allow for a good numerical 
illustration. Fig. [2] shows the evolution of the derivative 
fR,{R) of the resulting f{R) model, for the above choice 
of initial conditions. 



3 lAU ~ 1.496 X 10" m. 

* Numerical solutions in this paper are obtained using Mathemat- 
ica's differential and algebraic solvers, making use of tfie stiffness 
option as well as increasing the maximum step number when ap- 
propriate. Plots are also produced with the same software. 



At{R) 
dR 




FIG. 2: The derivative of the f{R) model for the initial 
conditions (37 1 (in Planck units). For large values of R, the 
model effectively behaves as B? gravity, since fniR) ~ R- For 
smaller R, the evolution enters the classical regime near the 
Gaussian fixed point. The Einstein-Hilbert term dominates, 
/rXR) becomes nearly constant, with a small positive slope 
reflecting the positivity of the scalaron mass-squared. At very 
small R (not shown) there is a departure from the Einstein- 
Hilbert action due to the IR divergence in the beta functions. 
This part of the action is never encountered as the system 
freezes at the IR de Sitter point. 



V. COSMOLOGICAL DYNAMICS 

Now, we proceed with studying the cosmological dy- 
namics of the model, i.e the cosmological fixed points, 
their stability and the transition from one cosmological 
era to the other, as well as if inflation can be viable in 
this scenario. 

A viable cosmological model, aiming to describe the 
background evolution of the universe from early to late 
times, should have a period of accelerated expansion at 
early times (inflation), followed by a radiation and mat- 
ter era respectively, evolving asymptotically towards de 
Sitter at late times. Each particular period has its own 
requirements in order to be viable. For example, a UV de 
Sitter point should be unstable, while an IR one stable, 
while the matter point should be a saddle with damped 
oscillation, so that strucutre formation has enough time 
to take place. 



A. Transition to the Einstein frame 

It will be useful for the latter analysis to first calculate 
the Einstein frame action, as an aid in calculating infla- 
tionary quantities, like the slow roll parameters. To do 
this we will introduce an auxiliary field and then confor- 
mally transform the metric appropriately. However, in 



the context of action ( 10 1 the latter transformation re- 



quires care, since Newton's G is running with curvature. 
Let us see this in more detail. Our original action 



(10) is a function of R,g{R),X{R) and implicitly of the 
metric through the Ricci scalar R. In the standard way. 



we introduce auxiliary scalars a and ^, and write our 
original theory in the Jordan frame as 



S 



1 

1 
2^ 



V^[ri<7)R-if'{a)<j^f{a))] (38) 
^[(t>R~V{^)], (39) 



with y(0) = (l)a{(j)) — f{a{(j)j), and (j) — /'(f)- We require 
that f"{cr) ^ 0, so that the function /' can be inverted 
to find cr as a function of (j). Note that the equation of 
motion for a gives the constraint which reproduces the 
original action, i.e a = R. 

The Jordan frame scalar </>, plays the role of the inverse 
of Newton's constant in front of R. For the transition to 
the Einstein frame, Newton's constant will have to be 
re-introduced through the conformal redefinition of the 
metric, and the question that arises in our scenario, is 



which Newton's constant should that be, since using a 
running G = G{R) could lead to ambiguities. We can 
resolve this issue by using Newton's G today, denoted 

We can now perform the conformal redefinition of the 
metric as 



9aP 



SttGo^ 



-9al3, 



combined with a redefinition on the scalar 



t>o exp 



'167rGo 



$ 



(40) 



(41) 



with 00 constant. Performing above two field redefini- 
tions in action ( 39 1 , we finally end up with the Einstein 
frame action 



S 



d^x 
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1 



V 167rGo 



R 



1 



(V$)^-L/($) +^„(g,7^,<i>). 



(42) 



The scalaron potential in the Einstein frame is then 
given in parametric form as, 



U{R) 



RfnjR) - fjR) 

^2 f^(U\2 ' V*"^'' 



2 (8^Go)' fR{Rf 



HR) = 



167rGo 



ln/fl,(i?). 



(44) 



The mass of the scalar $ is defined in the usual way 
through the Einstein frame potential as 



'■cff 






(45) 



We will use the above two relations later when we will 
work out the inflationary slow roll parameters. The Ein- 
stein frame potential for different values of p is plotted 
in Fig. [3J The maximum corresponds to the unstable de 
Sitter point in the UV, with the cosmological evolution 
occurring "to the left" of it, i.c to smaller field values. 
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FIG. 3: The Einstein frame scalar potential (in Planck 
units), described by relations (43l-(|44[|, for p = 0.8 (red, 
dashed), p = 1 (back, continuous) and p — 1.15 (blue, dotted) 
respectively. Cosmological evolution starts from the maxi- 
mum of the potential, which corresponds to the unstable UV 
de Sitter point, and evolves towards smaller values of the field 



B. de Sitter solutions 

Let us now look for the simplest cosmological solutions, 
which are the maximally symmetric constant curvature 
ones. In f{R) gravity, they correspond to the points 
where the potential in (14) has an extremum, i.e solution 



of the algebraic equation 

RJr{R) - 2f{R) = 0. 



One can check that the same condition is also derived in 
the Einstein frame by requiring that dU/d^ — 0. 
Using relation ( 20 ) , condition (|46| implies that 



2p^gpx + (1 - 2pA)/3g 



0. 



(47) 



For Px, Pg non-zero, and a given p, equation (47 1 de- 
fines a family of solutions, described by a curve in the 
g — A plane, which is the locus of all de Sitter points. 
(46) Any intersection of it with the RG trajectory will imply 



a de Sitter era in the particular cosmological evolution. 
It is interesting to note that any RG fixed point will al- 
ways satisfy the de Sitter condition (46) or (47), since 
there, (3g = P\ = 0. This is an identity for fixed points, 
as f{R) oc R^ there, but we can check that they are de 
Sitter by inspecting the Einstein frame potential. In par- 
ticular, the UV RG fixed point is a always a de Sitter 
point, as the potential (43 1 stays finite as R 



oo. 



The location of the de Sitter line depends on the value 
of the parameter p, which shifts the scale of both early 
and late time de Sitter points. As a starting point, we can 
get an idea of the de Sitter points structure by setting 
p = 1 in equation (47) and working out the resulting 
de Sitter line, which is shown in Fig. [ij The de Sitter 
line passes through the UV RG fixed point, yielding this 
way an infinite number of de Sitter points. This can be 
seen as follows: as pointed out before, the RG UV fixed 
point is a de Sitter point itself. On the same time, the 
behavior of the RG evolution in the vicinity of the RG 
UV fixed point is described by an unstable spiral, which 
circles the fixed point infinitely many times as fc — >■ oo 
(or R — >■ oo). As a consequence there will be an infinite 
number of intersections between the de Sitter line and the 
RG phase curve. The UV RG fixed point is the limiting 
de Sitter point of the above infinite set of de Sitter points. 



Furthermore, as Fig. [T] shows, there is an "outer" de 
Sitter point in the UV regime, and another one in the 
IR. For the case p = 1, wc find 



Ads)i 



WSj'^dSjUV 



(0.02,0.27), 



while it is easy to show that 



(gdS,Ads)iR- (0,0.25). 



(48) 



(49) 



Notice that the "inner" UV de Sitter points cannot be ac- 
cessed, since they are protected by the outer one. At least 
in the Einstein-Hilbert truncation, and under the cut-off 
identification considered here (fc^ = pR)i this seems to be 
a general behavior: There is always an infinite set of UV 
de Sitter points, all hidden by the most outer one, and a 
de Sitter point in the IR. As a consequence, classical cos- 
mological evolution cannot reach the extreme UV regime 
around the UV RG fixed point, i.e for k"^ = pR -^ oo. 

We now want to understand how the de Sitter line 
changes as we vary the dimensionless parameter p in our 
cut-off identification. There are two extreme cases lead- 
ing to two limiting de Sitter lines, one for p — >■ oo and 
another for p — >■ 0. Solving equation (47) for g and taking 
the limit p — > oo the limiting curve is described by 



5p^oo(A) = — f 12A2 - 4A - 3 + ^UAX-^ + 672A3 - 824A2 + 216A + 9 



(50) 



which for a realistic RG evolution gives a de Sitter point 
at the UV, and another one very close to the GFP, i.e 
X r^ g r^ 10""^^. Therefore, by tuning the parameter p 
to very large values, both UV and IR de Sitter points 
are shifted towards the UV along the RG trajectories. 
Notice that letting A —> in (50) we get gp^oo -^ 0, i.e 
the curve passes through the GFP at (A, g) = (0, 0). 
On the other hand, as p — > 0, the de Sitter line becomes 

5p^o(A) = ^(l-2A)2, 

which again gives a de Sitter point in the UV and a second 
one for (g. A) ~ (0,0.5), as g <C 1 in the IR regime, but 
now both points are shifted towards the IR. The de Sitter 
lines corresponding to the extreme cases described above 
can be seen in Fig. |4] 

To summarise: the general trend is that by making p 
smaller, the position of the UV de Sitter point is shifted 
towards smaller values of R (i.e moving away from the 
RG UV fixed point), while the situation is the opposite 
for increasing p. 



Let us turn attention to the stability of the de Sitter 
points. As said before, a de Sitter point is (un)stable if 




FIG. 4: The limiting de Sitter (continuous) and slow roll 
lines (dashed) respectively. Red colour corresponds to p — )• oo 
("bell" shaped curves), and the green to p — >■ respectively. 
For p — >■ CX3 both de Sitter and slow roll line go to zero as 
A^ 0. 



{m^g < 0) m^ff > 0. Therefore, the equation 
'm'lg{R,g,\) =0 

will in turn define a line on the g — X plane along which 
the square of the mass becomes zero. Another useful line 
on the g — X plane is the one along which the square 
mass diverges, i.e its denominator becoming zero. Both 



m 



off 



0, and m^g — > oo lines will divide the g — X plane 
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into regions of positive and negative mass squared. This 
can be seen in Fig. [T] for p = 1. If the UV de Sitter point 
should play the role of an inflationary era, it should be 
an unstable, while the IR one should be stable. We will 
see later that this can be achieved for a range of values 
for the parameter p. 

The general expression for the scalaron mass in terms 
of the dimensionless couplings is quite complicated, but 
it simplifies reasonably in the classical and IR regime, 
where g ^ 1- In this case, the mass takes the form 



off(IR) 



i?(2A-l)3 



365 (8AV + 4A2p + A(8 -6p) + p-2) 



, (51) 



where we neglected terms of order g^ and higher. The 
critical points where the denominator of the scalaron 
mass vanishes will signal a singularity, with the scalaron 
mass going to infinity. 

Let us study the positivity of Ti^ffdR-) by first studying 
the special case of the GFP regime, where in addition 
to g <C 1, it is also A ^ 1. In this case, relation (51) 
simplifies to 



R 



cff(GFP)-3g^(2-p)' 



(52) 



and the sign of it is positive when p < 2, while for p — !• 2 
it blows up. We recall that the renormalisation condition 
(13) fixes pRo/go ^ rrip, so the scalaron has a Planck- 
sca e mass near the Gaussian fixed point. 

In the IR regime, where A ~ 0(1), we have to study the 
full relation (51). We distinguish two regimes, one when 



p < 2 and another when p > 2. For p < 2, the vanishing 
of the denominator of (51 ) has only one relevant solution 



\ = \^ (p) being a function of p. 



< A < A* : mlg > 0, 
A* < A < 0.5 : mlg < 0. 



(53) 
(54) 



In the limiting case of p ^- 0, A* ^- 0, while as p ^ 2^, 
it is A* ~ 0.31. 

On the other hand, for p > 2, there are two relevant 
solutions, A^(i) and A*(2). We have the following cases 



< A < A*(i) : miff < 0, 
A*(i) < A < A^(2) : m^ff > 0, 
A*(2) < A < 0.5 : miff < 0, 



(55) 
(56) 
(57) 



with both A*(i), A*(2) varying with p, i.e A*(i) = A*(i)(p), 
A»(2)(p). In particular, we have that as p — > 2+ 



(2) 
*(1) 



and A»(2) — 0.31, while for p — !• oo A»(i) ~ 0.2 
and A*(2) — 0.5. 

We conclude from the above analysis that the case p > 
2 is rejected, since miff is negative around the GFP. On 
the other hand, for p < 2 the mass m^fj is positive around 
the GFP (A ^ 1) and it stays positive for A < A* with 
A* approaching A* ~ 0.3 as p — > 2. 

An important point when p < 2, concerns the position 
of the IR de Sitter point. From above, it is understood 



that both the position of de Sitter points as well as the 
critical point A*(p), beyond which m,eS becomes negative, 
depend on the parameter p. What it turns out to be is 
that the corresponding position of the IR de Sitter point, 
will lie ahead of A* on the A-axis for p < 0.9, which 
means that the de Sitter point will be unstable. As a 
result, all trajectories with .g ^ 1 and p < 0.9, will posses 
an unstable IR de Sitter point. In other words, the RG 
trajectory will pass through the mass singularity point 
A*, making m^fj negative, before the trajectory reaches 
its actual terminating (de Sitter) point. 

From the above stability analysis, we see that the pa- 
rameter p has been constrained to be 0.9 < p < 2. In 
the next section, we will further constraint p by requir- 
ing that the different cosmological periods are connected 
with each other in a viable way, finding that p ~ 1. 



C. Dynamical evolution from UV to IR 

We saw that in principle we can have de Sitter solu- 
tions, and the existence of a classical regime ensures for a 
standard radiation/matter era respectively. It is impor- 
tant though, that the cosmological eras are connected 
dynamically in a viable way. This will be the subject of 
this section. More precisely, we will consider action ( 10 ) 



in the presence of a perfect fluid with barotropic index 
w = p/p, and study its dynamics by means of a dynami- 
cal system analysis, by improving the dynamical system 
for f{R) gravity, presented in Ref. [5T], to account for 
our RG-inspired f{R) model. 

We can start by defining the following dimensionless 
variables. 



Xi 



X2 



X3 



X4 



HJr' 

-I 
6H^fR 

R 
6^' 

R^Pr 



(58) 
(59) 
(60) 
(61) 



with an over dot denoting derivative with respect to cos- 
mic time. The Hubble parameter is defined as H = d/a, 
with a the Universe scale factor. In the absence of radi- 
ation it is X4 — 0. 

Then, the background dynamics can be expressed in 
terms of the dynamical system |51| , 



x'f = -\- X2- 3X2 ^x\- XiX:i + X4, 


(62) 


A^ a;2(2a;3 - xi - 4), 


(63) 


4 = Z^i^_ 2x3(^3-2), 

TO 


(64) 


x\ = -2X3X4 + XiXA, 


(65) 
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with the constraint 



n,, 



~2 



1 — a;i — a;2 ~ 2^3 — X4, (66) 



and primes here denoting differentiation with respect to 
In a. 

The quantity m = m(r) is defined as 



dlnfR _ RJrr 
dlnR ~ fR ' 
d\nf _ RJr _ a^ 
d In i? / X2 



(67) 
(68) 



TO = m{r) characterizes the particular f{R) model, and 
it needs to be given a priori in order for the dynamical 
system to close. In principle, given a particular f{R) 
model, one is able to invert r — r{R) and plug into m 
to get m = m{r). However, in our case the form of the 
/(i?) model is dictated through the particular running of 
the couplings g{R), X{R), by solving the system of beta 
functions. Therefore, in order to close the dynamical 



system ( 62 )-( 65 ) we will need to evolve the couplings with 



time as well. 

In addition, the effective equation of state is given by. 



Woff = -^ (2x3 - 1). 
For the dimensionless couplings we can write. 



(69) 



^ ^ dRdrdN^ 2R dr Ux^^ + ^xg^ ' ' ^ ' 



, dXdR dr _ /3a dR ( dr , dr , 



A' = 



dRdrdN 2Rdr \dx2 fes ^ 



(71) 



After some algebra, we get 



P 



A' = 



where 



2R\flR-fRj-fRRfR 
/3a ( f 



2-^ V/r^ ~ /fl/ ~ JrrIR 



X3X2 ~ X^X-i 



X^X2 ~ X2XZ 



(72) 



(73) 



= x[{xi,g^\) through the relevant evolution 
equation. The complete dynamical set of equations is 
now (l62l)-(|65| supplemented with (fTSl-fTal). Notice that 



any fixed point of (62|-(65) automatically satisfies (72) 
( 73 1 as well. 



One should be reminded here that the derivatives 
with respect to i?, e.g Jr, can be explicitly expressed 
using (lw|, ([20]) and (ifl]). In addition, both r and 



TO are implicit functions of curvature i?, through r = 
r(A(7?),5(i?))and5 = 5(A(i?),5(i?)). 

The RG improved dynamical system with 0:4 = has 
three cosmological fixed points: An early time de Sitter, 
a matter, and a late time de Sitter point respectively. Of 
course, we expect that a radiation fixed point will appear 



by the time we introduce 2:4. For a complete analysis and 
the fixed point structure and their stability one can refer 
to Ref EI]. 

The de Sitter point Pi, the matter point P5 and the 
radiation point P7, are given in the general form P = 
{xi,X2,xz,Xi) as. 



Pi = (0,1,-2,0), 

, 3too 4too + 1 4too + 1 

Pf 



,0 



(74) 
(75) 



^ 'too + 1' 2(too + 1)2'2(too + 1) 

4too 2to,o 2?7io —brriQ — 2to,o + 1 



Too + 1 ' (toq + 1)2 ' Too + 1 ' (toq + 1)2 

(76) 

The de Sitter point Pi is characterised by r = —2, and 
is stable as long as 



< toL^„2 < 1- 



(77) 



On the other hand, the points P5 and Pg define a family of 
fixed points parametrized by to, all lying on the line to, — 
—r — 1. An acceptable matter era requires that standard 
GR is recovered, i.e to -> {/rr — 0), yielding P5 — 
(0, —1/2, 1/2), and therefore r = —1. For tti ~ 0, and in 
the presence of radiation, a radiation fixed point will also 
exist in the vicinity of P5. In particular, the existence of 
a saddle matter era requires that at the matter point, 



lr=-2 



+0, 



dm(r) 



dr 



> -1. 



(78) 



The shape of the curve to = m(r) on the m — r plane 
can provide us with sufficient information regarding the 
asymptotic behavior of the particular f{R) model. In our 
case, we can work out the m = m{r) curve by integrating 
the system of beta functions, and then evaluating both 
r = r{X,g) and g — g{\,g). By choosing a typical RG 
trajectory for p = 1 (i.e k^ = R), and initial conditions 



for the system of beta functions those of ( 37 1 , we get the 
m ~ r curve shown in Fig. [5] We see that cosmologi- 
cal evolution begins from an unstable (r > 1) early time 
de Sitter point, and then evolves towards the (radiation) 
matter point at (r, to) ~ (—1, 0). It then leaves the mat- 
ter point and evolves towards a stable IR de Sitter point 
at r = —2. Notice that the matter point is approached 



from positive values of m as condition ( 78 1 requires. 



For illustrative purposes. Fig. |6] shows the cosmological 
evolution from the matter to the IR de Sitter point in 
the coordinate space, while Fig. [8] shows the evolution of 
the effective index and slow roll parameter WeS and ey 
respectively, from the UV de Sitter point to the matter 
one. 
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FIG. 5: The m — r plane for p = 1 and the set of initial 



cond itions (37 1, with m{r) and r given by relations (671 and 
mm respectively. Point A corresponds to the unstable UV 
de Sitter point, point B to the saddle matter point, while C 
to the stable IR de Sitter respectively, as described in section 
|VC[ The dashed lines correspond to r = —2 and m — — r — 1 
respectively. 




FIG. 6: The c osm olo gica l trajectory described by the dy- 
namical system (621 - (651, in the space of the coordinates 



(xi, 12, sa), leaving the matter point and evolving towards the 
IR de Sitter. In particular, it spirals around the unstable mat- 
ter point, and then evolves towards the stable de Sitter in the 
IR. The initial conditions chosen are {Xo,go) ~ (10~^,10~^), 
and (a;io, 2:20,2:30) = (a;io(m) + 10"^,a:20(m) - 10~^,3:3o(m) + 
10~®), with (a:io(m), 2;20(m),a;30(m)) denoting the coordinates 
of the matter fixed point given i n (|75[ ), and mo is evaluated as 
7no(Ao,i7o) using ( |67| l and ( |20[ |-(|21p^ Above initial conditions 
give ro ~ —1.02, mo — 1.24 x 10~*. We also assumed that 
X4 — 0. The amplitude of the oscillation along the xi axis is 
of the order 10"^. 



As it turns out, under a suitable choice of initial con- 
ditions for g, A and p, it is possible to get a cosmology 
where the UV regime is correctly connected with the IR 
one. The question that arises is if there are any bounds 
on the parameter p in this direction. In fact, for p > 1.1 
the behavior of the evolution on the m — r plane starts 
becoming unstable, and evolution does not reach the late 
time de Sitter point, after leaving the matter era. What 
is more, as p increases the matter era happens to be ap- 
proached from negative values of ?ti, which as explained 
before is forbidden. Furthermore, as was also explained 
in the previous section, the positivity of ?7i^g in the IR 
regime (stability of IR de Sitter point) as well as in the 




Log(a) 



FIG. 7: The effective index w^s and the slow roll pa- 
rameter e (relations ( 69 1 and ( 82 1 respectively) from the 
UV de Sitter to matter domination for initial conditions: 
(a:io, 3:20, X'io, Ao, (?o) = (10-^ -1 - 10-^ 2 - 10-^ 0.26, 0.02) 
and In(ao) = —30, and e re-expressed as e = 2 — 3:3. 



GFP regime puts the extra restriction 0.9 < p < 2. 

Therefore, we conclude that the viability of both the 
classical regime and late time cosmology restricts p to lie 
in the range 

0.9</5<l.l. (79) 



D. Inflationary dynamics 

We showed that our particular f{R) model exhibits an 
unstable UV de Sitter point, which can be dynamically 
connected with the radiation/matter era in a viable way. 
We would like to understand if the UV de Sitter point, 
describing a primordial inflationary era, could be obser- 
vationally viable i.e if the scalar and gravitational fluctu- 
ations amplitudes as well as the number of e-foldings are 
those that are required according to observations. Recall 
that the only free parameter in our model is the dimen- 
sionless parameter p. 

Below, we will evaluate all inflationary quantities in 
the Einstein frame, ignoring the non-minimal coupling 
between matter and the scalar field, since inflation is a 
(almost) vacuum dominated period. 

Let us flrst revise some standard notions of scalar fleld 
inflation. To start with, the slow roll parameters ensure 
that the scalar field (inflaton) has a small kinetic energy 
during inflation, compared to the potential energy, so 
that the latter dominates. The two slow roll parameters 
are deflned as 



7? 



H 



d\nH 

dN ' 
1 de 



(80) 
(81) 



with the overdot denoting differentiation with respect to 
cosmic time. For a scalar field action with a kinetic term 
and a potential, they can be alternatively (and equiva- 
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lently to first order in e, 77) defined as 



ey 



W 



p 
IGtt 



IT 






(82) 
(83) 



with the subscript <& denoting differentiation with respect 
to the Einstein frame scalar field $ respectively. 

Inflation occurs as long as the slow roll condition is 
satisfied, i.e 



ev <C l,?7y "^ 1) 



(84) 



and ends when ev,ilv ^ 0{1). SniaUness of ey ensures 
that the spacetime during inflation remains sufficiently 
close to de Sitter, while smallness of r]v ensures that vari- 
ation of ev per e-fold is sufficiently small. 
The number of e-folds is given by 



N 



In^ 



U 



d$, 



(85) 



with Oj, Of the scale factor at the start and end of inflation 
respectively, and the slow roll approximation used in the 
last approximation. Above integral can be of course eval- 
uated in terms of the couplings and curvature R through. 



d$ = 



9$ dg 
~dg~dR 



9$ d\ 
~dXdR 



9$ 
dR 



dR, 



(86) 



and the integral (85) can be calculated between two 
points i?i and Rf along the RG trajectory. Notice that 
in the vicinity of a de Sitter point the number of e-folds 
diverges since there C/$ — >■ 0. 

Fluctuations of the scalar field during inflation, gener- 
ate scalar and gravitational perturbations, whose power 
spectra in the slow roll approximation are given by (see 

e.g. m) 



r,^ 



1287r C/3 



m^Ul 



(87) 



k=aH 



Va 



128 U 



iH 



assuming evaluation at the horizon crossing of the rele- 
vant mode. The scalar power spectrum becomes infinite 
when evaluated on a de Sitter point, reflecting the stan- 
dard infra-red divergence. This behavior can be seen in 

Fig.m 

Notice that expressing the derivative of the potential 
as 



dU 
9$ 



dU dR 
dRd^ 



16TTGa fa dU 
3 IrrBR' 



(89) 



and using a similar expression for the second derivative, 
the slow roll parameters in the Einstein frame can be 




FIG. 8: Upper row: The scalar (left) and gravitational 
(right) fluctuation power spectrum, as given by relations ( 92 1 
and (93 1 respectively, as a function of the couplings \,g, and 
setting p = 1. The scalar power spectrum peaks along the 
de Sitter line, as on a de Sitter point it is Ps — > 00. Lower 
row: The corresponding contour plots of the scalar (left) and 
gravitational (right) spectrum of upper row. In the scalar 
power spectrum the dotted line corresponds to the de Sitter 
line, along which the power spectrum diverges. Higher values 
correspond to lighter shaded areas. 



calculated to be 



ey(i?) = 



1 / 2f-RfR 
S\f-RfR 



wiR) = 3 



'^ fR + IrrIrR 



^/rrJ 



fRRiRlR - f) 



(90) 

(91) 



These relations can also be viewed as a function of the 
scalar $ 



$(i?), through relation (44), as well as func- 
tions of 5, A, p through relations (20) and (21). 



The equation ey(5, A) = 1 deflnes a curve in the g — 
A plane ("slow roll line"), whose intersection with the 
RG phase curve corresponds to the end of inflation, and 
is associated with the corresponding de Sitter line for a 
given p. The slow roll line for p — 1 can be seen in Fig. 

HI 

In general, decreasing p, the slow roll lines shift away 
from the UV RG fixed point along the RG evolution, 
and vice versa as p -^ 00. The opposite is true for the de 
Sitter lines, which means that an increasing p increases 
(decreases) the scale where inflation starts (ends), and 
the opposite is true for decreasing p. It is interesting to 
note that for p — ^ 00, the low energy de Sitter point lies 
before the point where ey = 1. The limiting slow roll 
lines for p — >■ 0, 00 are shown in Fig. [4J 



Let us now move to the power spectra, given in (87) 



13 



and (88). In order to match the scalar fluctuation am- 
plitude according to the CMB observations |52] , we need 
■Ps ^ 2 X 10"^ and Vg < 0.27's. The precise value of 
the amplitudes depends on a set of values for (g, A, p) 
evaluated at the particular scale of interest. It will be 
useful first to give the explicit expressions of the spectra 
in terms of g, A and p, for the beta functions (l6l [7|. We 
find 



-Ps^ 



^9 = 



128 A{g,X,p)B{g,X,pf 
3p C{g,\,p)^Dig,\,p)^' 
128Aig,X,p)Big,X,p) 



P C{g,\pY 

with the additional definitions 



(92) 
(93) 



A{gA,p)=g{Ag-{l-2\f), 

B{g, X, p) = 96g^p + g {{-24X^ + 4A + 6) p - 6) 

- (1 - 2XfXp, 
C(5, A, p) = -192.g2p + 4g (3 {AX'' - l) p + 2) 

+ (1-2A)2 
D{g, A, p) = -192g^p + 4g ((12A2 ~ 4A - 3) p + 4) 

+ (l-2A)2(4Ap-l). (94) 



We arrived at relations (92|-(|93[), using relations (89) and 



(20l-(21| to re-express the spectra appropriately. Anal- 
gous (but more complicated) expressions can be derived 
for beta functions with other gauges and cut-off func- 
tions. 

We have seen in the previous sections that stability 
requirements of the classical regime (GFP regime) as well 
as of the late time cosmology require that 0.9 < p < 1.1. 
Therefore, the first thing to investigate is inflation can 
be observationally viable for p in this range. 

So, let us proceed by studying the case of p = 1. In 
this case, we also know the values of the couplings at 
which inflation starts and ends , Pstart = (ffstart, Astart) — 
(0.02,0.27) and Pond = (gend,Ae„d) ^ (0.02,0.22), with 
-Pstart corresponding to the UV de Sitter point, and Pond 
to the point where ey = 1 (see also Fig. [l]). For the 
connection with observations one is in principle interested 
at the value of the power spectra about 60 e-foldings 
before the end of inflation. Now, for p = 1, and as can 
also be seen in Fig. |VD[ between Pstart and Pcnd both 
power spectra are smooth, decreasing functions of g and 
A, acquiring their lowest value at Pcnd, 



0.067, r„ ~ 0.052. 



(95) 



One sees that the (lowest) values of the power spectra 



(95), are too large to agree with observations, yielding a 



non- viable inflationary period for p = 1. It is not difficult 
to check that this behavior is true for all values of p 
between 0.9 ^ p ^ 1.1. Therefore, a viable late time 
cosmology cannot be combined with a viable primordial 
infiation. 



Having seen that an observationally viable inflationary 
era is not in agreement with a viable late time cosmol- 
ogy, which requires p ^ 1, we ask the following question: 
could inflation be viable on its own for some parameter 
p, away from p ^ 1? Let us try to understand this by 
checking the behavior of the power spectra ( 92 ) and ( 93 ) 
for the extreme cases of p — >■ and p — >■ oo respectively. 
Assuming a (finite) value of g and A we find that 



liiRVs.Vg == oo, 

p— ^0 



(96) 



which is obviously unacceptable. 

On the other extreme, i.e when p — ^ cx), the power 
spectra go to zero. 



lim Vs,Vg 

p— >oo 



0, 



(97) 



which is potentially viable. For the scalar to tensor ratio 
we find that 






p— J-OO 



48 (485^ + g (-12A2 + 4A + 3) - (1 - 2X)^Xy 

(96g2 + g (-24A2 + 4A + 6) - (1 - 2A)2A)^ 

(98) 



Remembering that when p 3> 1, the end of inflation, 
which is described on the phase space of g — X by the slow 
roll line, is shifted towards smaller values of the couplings, 
as can also be seen in Fig. [4] Therefore, we can get an 
estimate of above ratio by assuming that the fluctuations 
are produced at a point in the linear regime of the RG 
evolution, where g ^ X <^ 1, yielding Vg/Vs ^ 0(1), 
which is observationally unacceptable. 

Before concluding this section, let us comment on an- 
other possibility of understanding inflation in this sce- 
nario, that is modeling it as E? inflation ^53j using the 
f{R) model found in (33) at large R: 



I{R) 



Go 



{R - 2Ao) + 6(2 - p)pR^ 



(99) 



Matching to the perturbation amplitude, R^ inflation can 
account for the observations if the coefficient of the R? 
term is of order 10^^ [2 . Hence we see in approximate 
way how tuning p to very large values suppresses the 
perturbations. However, this results in an unacceptable 
classical limit as well as a non- viable late time cosmology 
for the reasons explained in previous sections. 

To conclude this section, it turns out that primordial 
inflation in this scenario cannot agree with observations 
unless p is very large, in which case the mass of the 
scalaron diverges and becomes tachyonic in the subse- 
quent evolution. Hence the observed fluctuations must 
be generated at a later period of inflation, which requires 
that more degrees of freedom should be introduced in 
the action, like for example a scalar fleld. Scalar fleld in- 
flation in the asymptotic safety scenario, and with scale 
identiflcation in the equations of motion, has been con- 
sidered in [55. A more exotic possibility is that the ex- 
tra degrees of freedom produce a flxed point with a very 



14 



small fluctuations. One notes that for small g, the tensor 
power spectrum becomes 



scalaron, provided 



I285A 



(100) 



which is suggestive that a fixed point with small g\ 
could be viable. Note the appearance of the product 
g\ ^ GA, which is the expected scale of tensor fluctua- 
tions in Einstcin-Hilbert gravity in a de Sitter phase with 
cosmological constant A. 



VI. DISCUSSION AND CONCLUSIONS 

We studied the cosmology of an f{R) model generated 
by the RG improvement of the Einstein-Hilbert action. 
The transition to f{R) gravity was achieved by identify- 
ing the renormalisation group scale to be proportional to 
scalar curvature. 



pR, 



(101) 



in the non-perturbative bet a- functions calculated from 
the exact renormalisation group equation. 

We found that the resulting f{R) model has some re- 
markable properties. Firstly, it maintains the correct sign 
for the graviton and scalaron kinetic terms. Very close to 
a non-trivial RG fixed point it behaves like R^ gravity, 
which is scale invariant, while it reduces to GR in the 
vicinity of the Gaussian fixed point. At solar and galac- 
tic scales, the scalaron's mass is of the order of Planck 
mass, preventing observable departures from GR at these 
scales. On the other hand, in the vicinity of the UV 
RG fixed point, the scalaron mass vanishes, reflecting 
the scale invariance of the action in that regime. 

The cosmological solutions of the f{R) model are also 
interesting. It naturally exhibits an unstable UV de Sit- 
ter point which evolves to a stable one in the IR. The 
effective cosmological constants are exponentially sepa- 
rated when Newton's G and the cosmological constant 
are matched to their observed values. What is more, 
there are an infinite set of de Sitter points as the UV RG 
fixed point is approached {R — > 00). However, classical 
cosmological evolution starts from the outermost de Sit- 
ter point, and therefore the UV RG fixed point is hidden 
behind it, and cannot be accessed. The Big Bang sin- 
gularity is avoided, since the de Sitter point is reached 
at infinite time in the past, i.e as i — >■ — c». The model 
therefore satisfies the requirements of a successful f{R) 
model itemized in Ref. |57l. 

Introducing matter content to the cosmology, we found 
that the UV de Sitter point can be connected to the IR de 
Sitter era through a radiation/matter era, with a stable 



0.9 <p< 1.1. 



(102) 



Unfortunately, the fluctuations generated during infla- 
tion at the outer UV de Sitter poi nt are too large to 
account for the observations (Section VD|. 

Therefore observable inflation requires extra degrees 
of freedom in the action, for example a scalar field driv- 
ing inflation at a lower scale. A more remote possibility 
would be that the extra degrees of freedom move the flxed 
point to a smaller value of gX, which could suppress the 
fluctuations. 

To make the comparison with previous cut-off identi- 
fications in the literature, performed at the level of the 
equations of motion, our constraint for the parameter p, 
i.e /9 '^ 1, is broadly consistent with scale identifications 
made in the equations of motion, rather than the action 
as here. In particular, in Ref. |40j . it was numerically 
found that for the identification k'^ ~ cH^, the constant 
c should be of order one, which is consistent with p ^ 1 
in our identification. 

This model can be improved by extending the analy- 
sis performed in this paper to higher truncations, i.e by 
including higher order curvature terms in the action. It 
is interesting to ask what features are generic. The exis- 
tence of a UV fixed point seems to be a universal feature 
of all truncations found so far, so we expect the Einstein 
frame potential of the scalaron to tend to a constant at 
large values of the field. However, we do not expect the 
presence of an inflnite number of de Sitter points to be 
generic, as it arose from the complex eigenvalues of the 
fixed point, which are not present for the general four- 
derivative truncation |55| . We should also include matter 
fields in the renormalisation group equations. With these 
modifications it might turn out that there is a model for 
which both early and late time cosmology agrees with 
observations. 

Note added. While this work was being finalised, Ref. 
[55] appeared, which also makes the scale identification 
in the action, and flnds an infinite number of de Sitter 
solutions near the UV fixed point. 
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